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Abstract 

We study, both analytically and numerically, the phenomenon of energy dissipation in single¬ 
domain ferromagnetic nanoparticles driven by an alternating magnetic field. Our interest is focused 
on the power loss resulting from the Landau-Lifshitz-Gilbert equation, which describes the preces- 
sional motion of the nanoparticle magnetic moment. We determine the power loss as a function of 
the field amplitude and frequency and analyze its dependence on different regimes of forced pre¬ 
cession induced by circularly and linearly polarized magnetic fields. The conditions to maximize 
the nanoparticle heating are also analyzed. 

PACS numbers: 75.50.Tt, 76.20.+q, 84.60.-h 
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I. INTRODUCTION 


Single-domain ferromagnetic nanoparticles are of great interest due to their unique 
physical properties such as superparamagnetism,— ^ macroscopic quantum tunneling of 


magnetization, 


3,4 


size 


-dependent characteristics,- and exchange bias.- 6 ' 1 These and other 


nanoscale properties of ferromagnetic nanoparticles make them very attractive for ap¬ 
plications, e.g., in high density data storage,- - — spintronic devices,— - — and biomedical 


engineering. 


14-17 


If the magnetic state of nanoparticles is controlled by a time-dependent 


external magnetic field, then nanoparticles absorb energy from the field and heat up. While 


the heating of nanoparticles is undesirable for most applications, this proper 


importance for magnetic hyperthermia applications- 


14 . 16.17 


'see also Refs. 


18 


y is of crucial 


20 and refer¬ 


ences therein). In a ferrofluid subjected to an external periodic magnetic field two heating 
mechanisms are usually considered,— 1 ^ one of which is related to Brownian rotation of 
nanoparticles and the other to their Neel relaxation. With these thermal-induced mecha¬ 
nisms, the energy dissipation per period is expressed in terms of the imaginary part of the 
magnetic susceptibility of the ferrofluid. It should be noted that since the magnetic sus¬ 
ceptibility is a function of only the Brownian and Neel relaxation times,— the deterministic 
dynamics of the magnetic moment of nanoparticles does not affect the energy dissipation. 


In contrast, if the rotation of nanoparticles is not allowed and the superparamagnetic state 
is not realized (this occurs, e.g., when the nanoparticles are embedded in a solid matrix and 
the temperature is small enough), then the heating phenomenon is expected to be strongly 
dependent on the dynamics of the nanoparticle magnetic moment. On a phenomenological 
level, it can be described by the deterministic Landau-Lifshitz-Gilbert (LLG) equation.— 1 ^ 
Because of its nonlinearity, the precessional motion of the magnetic moment can be very 
complex. In particular, the circularly polarized magnetic field, whose polarization plane is 
perpendicular to the anisotropy axis, can generate the periodic and quasiperiodic regimes of 
precession of the magnetic moment.— - — Moreover, the precessional motion of the magnetic 
moment induced by the linearly polarized magnet field can exhibit chaotic behavior.— - — 

The aim of this paper is to study the dependence of energy dissipation on these regimes 
of precession of the nanoparticle magnetic moment. To the best of our knowledge, this 
problem has not been addressed before. We emphasize that the above-mentioned regimes 
and transitions between them can exist only in anisotropic nanoparticles. Note that the 
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energy dissipation in such nanoparticles, arising from the precessional motion of the magnetic 


moment, was the subject of Ref. 


32 


But the authors considered only the periodic regime 


of precession induced by the circularly polarized magnetic field, ft is also important to 
stress that, because the influence of eddy currents on the magnetic moment dynamics can 
be accounted for by introducing an additional damping parameter,— the analysis presented 
below is applicable to both dielectric and metallic nanoparticles. 

The paper is organized as follows. In Sec. HU we describe the model and introduce the 
reduced power loss. The analytical solutions of the LLG equation, obtained in the small 
amplitude approximation for both circularly and linearly polarized magnetic fields, and the 
corresponding power losses are presented in Sec. IIII1 Our numerical results are reported in 
Sec. mu Here, the reduced power loss and its connection with the character of the precessional 
motion of the nanoparticle magnetic moment are studied depending on the amplitude and 
frequency of these magnetic fields. We summarize our findings in Sec. |Vj 


II. DESCRIPTION OF THE MODEL 

We consider the Stoner-Wohlfarth particle^ of spherical shape characterized by the uni¬ 
axial anisotropy held H a and the magnetic moment m = m(t) with a constant magnitude 
|m| = m. It is assumed that the z axis of the Cartesian coordinate system xyz is directed 
along the particle easy axis and the magnetic moment m is under the action of both the 
alternating magnetic held 


h(£) = hcos(ut)e x + phsm(u>t)e y (2.1) 

and the static magnetic held H = He z . Here, e^, e y , and e 2 are the unit vectors along 
the corresponding axes of the Cartesian coordinate system, h and c o are the alternating 
held amplitude and frequency, respectively, and p = —1, -|-1 or 0. The case with p — ±1 
corresponds to the circularly polarized magnetic held h(£) rotating in the xy plane in the 
clockwise (if p = —1) or counterclockwise (if p = +1) direction, and h(t) is linearly polarized 
along the x axis when p = 0. The magnetic energy of such magnetic moment is given by 

W = -^-m 2 — m • H — m • h(£) (2.2) 

2m 

(m u = m • e„, v — x, y, z) with the dot denoting the scalar product. 
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We first assume that the dynamics of the nanoparticle magnetic moment is governed by 
the stochastic LLG equation 2 ’?-- 

dm . a dm 

— = — 7 m x (H ef r + n) H-m x —. (2.3) 

at m at 

Here, y(> 0) is the gyromagnetic ratio, a(> 0) is the dimensionless damping parameter, 

r)\Al ty) 

H cft = = H a — e, + H + h (t) (2.4) 

0 m m 

is the effective magnetic held acting on the magnetic moment m, and the cross sign denotes 
the vector product. As usually, the Cartesian components n v (t) of the thermal noise n = n(i) 
are considered as independent Gaussian white noises characterized by zero means, ( n u {t)) = 
0, and correlation functions (n l/ (ti)n l/ (t 2 )) = 2A S{t- 2 — t\), where the noise intensity A is 
proportional to the thermal energy ksT (k^ is the Boltzmann constant, T is the absolute 
temperature), 5(t) is the Dirac 5 function, and the angular brackets denote averaging over 
all realizations of n v (t). I 11 general, due to the thermal fluctuations, the dynamics of m is 
stochastic, as in the case of nanoparticles in the superparamagnetic state. However, under 
certain conditions (see below) the nanoparticles can be single-domain and, at the same time, 
their magnetic moment dynamics can be approximately described by the deterministic LLG 
equation, i.e., Eq. (12. 3 ll with n = 0. Our purpose is to determine the power loss under these 
conditions. 

The condition that nanoparticles are in the single-domain state follows directly from the 
Brown’s fundamental theorem.— Since, according to it, the single-domain state is energeti¬ 
cally favorable if the nanoparticle diameter d is less than a critical value d max (which is of 
the order of the domain wall thickness), this condition can be written as d < d max . Next, 
the thermal fluctuations do not play an important role in the dynamics of m if the ther¬ 
mal energy fcgT is much smaller than the smallest energy scale in Eq. (12. 2 [) . Because the 
condition h = h/H a <C 1 is assumed to be realized, this energy scale is given by mH a h or 
2 KVh, where K = H a M /2 is the uniaxial anisotropy constant, M = m/V is the nanopar¬ 
ticle magnetization, and V = 7 rd 3 /6 is the nanoparticle volume. From this it follows that 
the thermal energy can be neglected if k — 2KVh/(kBT) 1. Introducing the parame¬ 

ter d\ = [3kBT/(rcK /r)] 1 / 3 , interpreted as the nanoparticle diameter for which k — 1 , we 
can rewrite the condition k 1 in the form (d/d\) 3 1. The last inequality is satis¬ 

fied with a good accuracy if d > d min , where d min can be chosen, e.g., as d min = 3d\ (in 
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this case min k = 27). Thus, the nanoparticles with d e (c Lin , d m ™) are single-domain 
and their magnetic dynamics is almost deterministic. Note that the interval (d m i n , d max ) 
exists (i.e., d min < d max ) if T < T max , where T max = 2KV max h/(27k B ) is the characteristic 
temperature, defined as the solution of the equation k\ d=d max = min k with respect to T, 
and V max = 7rd^ ax /6. In other words, for each single-domain nanoparticle there is a finite 
temperature interval (0,T max ) in which the thermal energy is negligible. 


It is important to stress that even if the condition k 1 holds, there always exist the 
thermal fluctuations of m leading to a significant (of the order of 2 KVh or greater) change 
of the magnetic energy. But if the average time interval (t) between these fluctuations 
essentially exceeds the calculation time (it can be chosen as 2ttN/u with IV 1), they do 
not influence the dynamics of m if oj cuo, where ooq = 2nN/ ( t ) is the characteristic thermal 
frequency. Associating (t) with the mean first-passage time for the magnetic moment,— (t) = 
\Jit/ne K (2auj r )~ l (ay = 'yHa is the resonance frequency), we obtain uy = 4i/7m e~ K au r N. 
So, if d e (dmin, d max ) and cu 3> oJo then the nanoparticles are single-domain and the dynamics 
of their magnetic moments can be considered as deterministic. We emphasize that these 
conditions are not too restrictive. In particular, according to Ref. |37|, the Co nanoparticles at 
room temperature T = 300 K are characterized by the parameters K = 4.12 x 10 6 erg/cm 3 , 
47 tM = 1.79 x 10 4 G, and d max = 96.4 nm. Therefore, assuming that h = 0.1, from the 
definition of d m ; n one obtains d min = 13.7nm. Then, choosing d = 15 nm, a = 0.1, N = 10 6 , 
and taking into account that in the considered case u> r = 10 11 s _1 , we find k = 35.2 and 
ujo = 2.2 x 10 2 s -1 . Since cuo strongly decreases with increasing d (e.g., for d = 17 nm we 
have k = 51.2 and ojq = 2.9 x 10 _5 s^ 1 ), there is almost no restriction on the alternating 
field frequency uj. 


Thus, the above estimations clearly show that if the single-domain nanoparticles are not 
too small then the description of the dynamics of the nanoparticle magnetic moments by 
the deterministic LLG equation is quite justified, and this approach can be used even at 
room temperatures. Since we restrict ourselves to this case, below the dynamics of m is 
considered as purely deterministic. 

In spherical coordinates, the deterministic LLG equation (12.31) (when n = 0) reduces to 
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a system of two ordinary differential equations 

(1 + a 2 )0=— a sin 6Jcos 0 + H) + ah cos OF + h F v , 

(1 + a 2 )tp =cos 0 + H — h cot OF + ah esc 0 TJ, (2.5) 


where 0 = 0(i) and cp — cp(t) are the polar and azimuthal angles of the vector m, respectively, 
t = u r t is the dimensionless time, the overdot denotes differentiation with respect to t, 
H = H/H a , 


F = cos p cos (cut) + psiiup sin(aA), 


( 2 . 6 ) 


u = u/u r , and F v = dF/dp. 

The dynamics of m is accompanied by the dissipation of magnetic energy W. The power 
loss, i.e., the magnetic energy dissipation per unit time, is defined as Q = lim r _ > . 00 (l/r) / Q r dtq, 
where q = —dW/dt is the instantaneous power loss. Since according to Eqs. (12.2ft and (j2.4[) 
q = H eff • dm/dt, the reduced power loss Q = Q/(H a mco r ), which is the quantity of our 
main interest, can be written in the form 


Q 


Jim —— 

t— too 7TIT 


dtn 


eff 


• Ill 


(2.7) 


with H e ff = H e ff/i7 a and r = co r r. To calculate Q , we need to solve the LLG equation that, 
in general, can be done numerically. But in some special cases the expression for the power 
loss can be determined analytically. 


III. POWER LOSS: ANALYTICAL RESULTS 

A. Circularly polarized magnetic field 

There are two qualitatively different regimes of the steady-state dynamics of m in the 
circularly polarized magnetic field rotating about the nanoparticle easy axis, namely, peri¬ 
odic and quasiperiodic.~ The analytical results are mainly available for the periodic regime 
characterized by the constant precession and lag angles, 0 = lim^^ 0 and <f> = lim^^ 0, 
where = p — put. As it follows from the system of equations (I2.5jl , the precession angle 
satisfies the equation^— 

.—- 1 — (“) 

h 2 = -———[(cos© + H — pu) 2 + (cub cos©) 2 ] (3.1) 

cos- © 
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and the lag angle is connected with the precession one as follows: 


. au . 

Sill <f> = -O^r- Sill 0. 

h 


(3.2) 


It has been shown (see also Ref. 


38]) that if the direction of held rotation is opposite to 


the direction of the natural precession of m then the periodic regime is always stable with 
respect to small perturbations. In contrast, if these directions coincide then in the parameters 
space there are the regions where a given periodic regime becomes unstable. Depending on 
the system parameters and initial conditions, the instability leads to the transition of the 
magnetic moment into one of three possible steady-state regimes. Two of them are periodic 
and one is quasiperiodic. In one of these two periodic regimes the sign of rn z is the same 
as in the given periodic regime and in the other it is opposite. The transition of m into 
the periodic regime with opposite sign of rn z corresponds to the irreversible switching of the 
magnetic moment (for more details see Sec. UVD . 

Integrating by parts, Eq. (12.7P can be represented in the form 


Q = lim i 

771 t— >oo T 


dt[m y cos (cut) — pm x sin(u;£)]. 


(3.3) 


From this, writing m x and rn y in spherical coordinates and using the relation p = <f> + put 
and Eq. (13.21) . we obtain the following general expression for the reduced power loss in the 
case of periodic regime: 

Q = au 2 sin 2 ©. (3.4) 


This quantity can easily be calculated at h -C 1. Indeed, introducing the designation 
a = sgn(cos@), where sgn (x) is the sign function, and assuming that 1 + aH > 0 (this 
assumption does not restrict the generality of the expression below), from Eq. (13.IK one gets 


0 


7T 


(l-a) + 


ah 


(1 + aH — apu) 2 + (au) 2 


(3.5) 


Therefore, since in the limit of small rotating held amplitudes the condition sin© = h[(l + 
aH — apu) 2 + ( au ) 2 ]^ 1 / 2 holds, Eq. (13.41) reduces to 


Q = 


a 


h 2 u 2 


(3.6) 


(1 + aH — apu) 2 + (au) 2 
According to this result, if ap = —1 or, in other words, if the direction of held rotation 
and the direction of the natural precession of m are opposite, then the reduced power loss 
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is a monotonically increasing function of u with Q ~ ah 2 uj 2 /( 1 + oH) 2 as u —> 0 and 
<5|w=oo = oth 2 /(l + a 2 ). In contrast, if op = 1, i.e. , if these directions coincide, then Q is a 
unimodal function of uj, which at u — o>o, where 

ujq = 1 + oH , ( 3 - 7 ) 

possesses an absolute maximum with 

h 2 

Q\co=u 0 — —• ( 3 - 8 ) 

a 

At the same time, the small-frequency behavior of Q and the limiting value Q |c=oo are 
the same as in the previous case. It should also be noted that the frequency ujq, at which 
the power loss reaches the maximum, is always larger than the resonance frequency o; res = 
(1 + oH)/(l + a 2 ), at which the precession angle ( 13 . 51 ) becomes maximal (if o = 1) or 
minimal (if o — — 1). 


B. Linearly polarized magnetic field 


If the alternating magnetic field (12. ip is linearly polarized (i.e., p = 0) and its reduced 
amplitude h is small enough, then the solution of Eq. (12. 3 p with n = 0 can be represented 
as m = ome z + mj_. In the linear approximation in h , we have m± z = 0, 

H e ff = (a + H)e z + hcos(ut)e x , (3.9) 


and thus the deterministic LLG equation (12.31) reduces to 

rii^ = (o + H)e z x m_|_ + aae~x riij_ — omh cos(u}t)e y . (3.10) 

Since in this case mi = mi^e^ + m± y e yi Eq. (13.101) is equivalent to the system of two 
equations 


= — (cr + H)m± y — oarhj_ y , 

(3.n) 

m i y = (cr + H)m± x + oam_ Lx — omh cos (ojt), 
whose steady-state solution can be found in the form 

m_ix = m[acos(cot) + 6sin(a3t)], 

m± y = m[ccos(ojt) + dsin(cbt)]. 
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(3.12) 








Substituting (j3. 12[) into (13. lip and using the linear independence of the trigonometric func¬ 
tions sin (cot) and cos(ujt), one straightforwardly obtains 

(1 + aH) [(1 + a H) 2 — uj 2 + a 2 uj 2 ] 


a = 


b = 


c = 


h, 

h, 


d = 


[(1 + aH) 2 — uj 2 + a 2 uj 2 ] 2 + 4 a 2 uj 4 
acu[(l + aH) 2 + uj 2 + a 2 uj 2 ] 

[(1 + aH) 2 — uj 2 + a 2 uj 2 ) 2 + 4 a 2 u) 4 
2era(l + aH)uj 2 

[(1 + aH) 2 — uj 2 + a 2 u> 2 ] 2 + 4a 2 uj 4 
cru;[(l + aH) 2 — uj 2 — a 2 uj 2 \ 


(3.13) 


-h. 


h. 


[(1 + aH) 2 — uj 2 + a 2 u 2 ] 2 + 4 a 2 u 4 
According to Eq. (j3.9|) . the reduced power loss (12.71) in the reference case becomes 


Q = — Jim h 

m t —>oo t 


dt cos(ujt)m± :L 


(3.14) 


From this, in accordance with Eqs. (I3.12p and (13. 131) . one gets Q = ujbh/2 and 

q _ ®h 2 uj 2 [( 1 + aH) 2 + u 2 + a 2 uj 2 } 

^ 2 [(1 + aH) 2 -u 2 + a 2 u 2 ] 2 + 4a 2 u 4 ' 

A simple analysis of the frequency dependence of Q shows that Q rs j ah 2 uj 2 /[2(1 + aH) 2 } as 
uj — >■ 0, Q\ti =00 = a/r 2 /[2(l + a 2 )}, and if a > \/3 then Q monotonically increases with uj. In 
contrast, if a < \/3 then Q is a nonmonotonic function of uj, which at uj = ujq, where 


UJq — 

reaches the maximum value 


3 — a 2 


1 + 


\/l + a 2 


(1 + aH), 


Q\u=L)q 


ah 2 


\/l + a 2 (2 + Vl + ol 2 Y 


(3.16) 


(3.17) 


4 (\/l + a 2 — 1 + a 2 ) 2 + a 2 (2 + \J 1 + ck 2 ) 2 
Comparing the power losses in circularly and linearly polarized magnetic fields, Eqs. 
and (I3.15p . we see that in the former case Q |s<i and Q\^ =00 are two times larger than in 
the latter one. Similarly, comparing Eqs. (13. 8 p and (I3.17p . one can make sure that Q\z=z 0 at 
a 1 is four times larger. The typical dependencies of Q on uj, calculated using Eqs. (I3.6p 
and (13.151) . are shown in Fig. [H The case with a < \/3 is illustrated by Fig. [Hr, and the case 
with a > y/3 by Fig. [ljo. As seen, the power loss in nanoparticles driven by the circularly 
polarized magnetic held exceeds that for the linearly polarized held of the same frequency 
and amplitude. It can be thus concluded that in the small amplitude approximation the 
nanoparticle heating is more efficient in the circularly polarized magnetic held. 
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FIG. 1. (Color online) The reduced power loss Q as a function of the reduced frequency u in the 
small amplitude approximation. The parameters used in Eqs. (13.61) and (13.15|) are h = 0.1, H = 0, 
a = +1, p = +1, a = 0.1 (a) and a = 2 (b). 


IV. POWER LOSS: NUMERICAL RESULTS 


Due to the nonlinearity of the LLG equation, with increasing the alternating field am¬ 
plitude the dynamics of the magnetic moment m can become very complex. In particular, 
the regimes of quasiperiodic— and chaotic^- 1 motion of m can be realized in circularly and 
linearly polarized magnetic fields, respectively. Therefore, to find the power loss in these 
and other cases, Eqs. (12.51) should be solved numerically in a wide region of the reduced 
amplitudes and frequencies. To this end, one needs to perform a number of runs for each 
value of h and u) from the intervals (h min , h max ) and (cu min , D max ) with some steps Ah and 
Auj. It is important to note that the transitions between different dynamical regimes of m 
can be irreversible and can depend on the trajectory in the space of discrete variables h 
and u, as it was shown for the circularly polarized field. 3910 Therefore, to avoid confusion, 
we first fix uj and then change h from h m i n to h max with the step Ah. One run consists 
in finding, for given h and uj, the solution of Eqs. (12.5p on the time interval (0,f sim ). We 
use the fourth-order Runge-Kutta method and consider the simulation time £ sim to be much 
larger than the time fi oss during which the magnetic moment, exhibiting regular motion, 
loses memory of its initial orientation defined by the angles 9o = 9(0) and cpo = In 

this case, we choose 9 0 = 10 -4 , cp 0 = 0, t sim = 3 • 10 4 , ti oss = 10 2 , and associate the reduced 
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power loss (12.7[) with the numerically obtained result 



dt[—(H + cos 6) sin 6 6 


+ hF cos 6 6 + hFy sin 6 p\. 


(4.1) 


In the chaotic regime, the solution of Eqs. (12.5ft is sensitive to initial conditions. There¬ 
fore, to avoid the dependence of Q on these conditions arising from the finiteness of t sim , 
an additional averaging of Q over (po is performed assuming that this angle is uniformly 
distributed in the interval (0,27 t). 

A. Small amplitude case 

Using the numerical procedure described above, we first analyze the dependence of the 
power loss (14. ip on the alternating held amplitude for rather small values of h. If the held 
frequency u is also small then Q as a function of h is well described by Eq. (13. 6 p (in the 
case of circularly polarized magnetic held) or Eq. (I3.15p (in the case of linearly polarized 
magnetic held). But if uj is relatively large then Q undergoes a qualitative change as h 
increases. Such behavior of the reduced power loss is demonstrated in Fig. [2] As seen, the 
analytical and numerical results are almost identical for small amplitudes and, at the same 
time, the difference between these results becomes very pronounced even for not too large 
values of h. The most remarkable feature of Q as a function of h is the existence of critical 
amplitudes, h r p! and h l £ for circularly and linearly polarized magnetic helds, respectively, 
at which Q changes abruptly. In the case of circularly polarized held, this corresponds to 
the so-called P-P transition from one periodic regime to another—— (see also Fig. |3|). The 
fact that Q\j 1= j 1 cp_ £ < Q\h = h'^ +£ ( £ 1) is a consequence of Eq. (13.4p and the condition 

@| j i= jcp_ e < 0|^c ?+£ . The similar transition occurs also in the linearly polarized held at 
h = h'p. . But in contrast to the previous case, the precession angle depends on time both 
before and after transition. 

Because the behavior of m in circularly and linearly polarized magnetic helds with arbi¬ 
trary values of uj and h is quite different, we consider these cases separately. 
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FIG. 2. (Color online) The reduced power loss Q as a function of the reduced amplitude h for 
circularly (p = +1) and linearly ( p = 0) polarized magnetic fields. The solid and dashed lines 
represent the theoretical results obtained from Eqs. (13.61) and (|3.15l) . respectively, and the symbols 
show the numerical results obtained from Eq. (14.11) . It is assumed that a = +1, a = 0.1 and 
u = 0 . 8 . 

B. Circularly polarized magnetic field 

Using Eq. (14. ip . we numerically calculated the reduced power loss Q for a wide region in 
the space of parameters uj and h. The results for this quantity and the boundaries between 
different periodic (P) and quasiperiodic (Q) regimes of the steady-state precession of the 
magnetic moment are shown in Fig. [3] for H = 0. Region 1 represents the periodic regime 
of precession of m, which is described by Eqs. (13. lj) and (13. 2 j) with 0 < 7 t/ 2 (i.e., a = +1) 
and p = +1. When h increases, this precession becomes unstable^ and, depending on 
the reduced frequency, the magnetic moment can make a transition to one of three steady 
states. These states are respectively characterized by (i) periodic precession with 0 < 7 t/ 2 
(region 2), (ii) periodic precession with 0 > 7 t/ 2 (region 3), and (iii) quasiperiodic precession 
(region 4). The precession angle 0 as a function of h is discontinuous at the boundaries 
(denoted by circles) between regions 1 and 2 (@|i < 0| 2 < 7 t/ 2) and between regions 1 and 
3 (@|i < 7 t/ 2 < 0| 3 ). An important difference between these P-P transitions is that the 
former is reversible, while the latter, corresponding to the magnetic moment switching, is 
irreversible. 

In region 4, the steady-state precession of m is quasiperiodic, i.e., the angles 6 and </> = 
(p — put are periodic functions with the same period, which in general is not commensurable 
with the field period 2ir/u. The transition from the periodic regime of precession in region 
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FIG. 3. (Color online) Color map for the reduced power loss Q and the diagram for the steady- 
state regimes of precession of the magnetic moment driven by the circularly polarized magnetic 
field. The regions with different periodic regimes of precession are indicated by numbers 1-3, and 
the region with the quasiperiodic regime of precession by number 4. The numerical results are 
obtained for p = +1, H = 0 and a = 0.1. 



FIG. 4. (Color online) Examples of steady-state trajectories of the magnetic moment driven by the 
circularly polarized magnetic field. The simulation parameters are chosen to be p = +1, H = 0, 
a = 0.1 and = 0.8. Trajectories for periodic regimes in regions 1 (h = 0.05), 2 (h = 0.1) and 
3 (h = 0.4), and for the quasiperiodic regime in region 4 (h = 0.26) are shown in (a) and (b), 
respectively. 


1 or 2 to the quasiperiodic regime in region 4 (P-Q transition, triangle line) is reversible. In 
contrast, the transition from the quasiperiodic regime of precession to the periodic regime 
in region 3 (Q-P transition, square line) is irreversible. As an illustration, in Fig. [4] we show 
the steady-state trajectories of m in regions 1-3 (a) and in region 4 (b). 

The numerical results for the reduced power loss Q are show in Fig. [3] as a color map. 
For the periodic regimes (regions 1-3), these results are in excellent agreement with those 
obtained from Eq. (13 . 4p . We note the following features of the reduced power loss. First, 
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FIG. 5. (Color online) The reduced power loss Q as a function of the reduced amplitude h of circu¬ 
larly polarized magnetic field for different values of the reduced magnetic field H. The parameters 
p, a and u are the same as in Fig. [4j 

since ©|i < ©| 2 , the P-P transition to region 2 is accompanied by an abrupt increase of Q. 
Second, the transition to region 3 is followed by an abrupt decrease of Q (see also inset in 
Fig. [3]). Because after transition to region 3 the direction of the natural precession of m 
becomes opposite to the direction of held rotation, this occurs for both the P-P and Q-P 
transitions. Third, the P-Q transition from region 1 or 2 to region 4 does not lead to a 
discontinuity in Q. And fourth, as is clearly seen from this figure, the maximum of the 
reduced power loss is reached near the (triangle) line of the P-Q transition. 

It should also be noted that the static magnetic held changes the magnetic moment 
dynamics and hence inhuences the power loss. For illustration, in Fig. [5] we show the 
dependence of Q on h for different values of H. The jumps of 0 at II = 0 correspond to the 
P-P (reversible) and Q-P (irreversible) transitions, while the jumps at H = ±0.5 correspond 
to the P-P (irreversible) transitions associated with switching of m. 


C. Linearly polarized magnetic field 

In this held, the dynamics of the magnetic moment differs considerably from that de¬ 
scribed in the previous section. One of the differences is the absence of the periodic regime 
of precession of m in the above sense. But the most striking difference is that the linearly 
polarized magnetic held can induce the chaotic regime of precession. 1 This implies the 
existence of regions in the parameter space in which the regular dynamics of m is still very 
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FIG. 6. (Color online) Examples of regular (a) and chaotic (b) trajectories of the magnetic moment 
driven by the linearly polarized magnetic field. The simulation parameters are p = 0, H = 0, 
a = 0.1, oj = 0.8, h = 0.4 (a) and h = 0.46 (b). 


complex (pre-chaotic behavior). For comparison, in Fig. [6] we plot the trajectories of m 
performing regular and chaotic (in a finite time interval) precessions. 

Using the previously described procedure, we can compute the reduced power loss Q 
in the case of linearly polarized magnetic held as well. To analyze its dependence on the 
character of precession, we need to establish whether the precession is regular or chaotic for 
a given set of parameters. This can be done by determining the sign of the largest Lyapunov 
exponent Ai (Ai > 0 corresponds to the chaotic behavior), which describes the divergence 
of neighboring trajectories.—^2 For our system, this quantity can be introduced as 

k 

Ai = lim \~ In y/\59(n At)] 2 + \5<p(n Ar)] 2 , (4.2) 

fc-s-oo k Ar ' 
n= 1 

where Ar is a short time interval (At <C Ar <C £ sim ), At is the simulation time step, and 56 
and 5(p satisfy the system of linear equations 


(1 + a 2 )59 = f e 59 + fpStp, 
(1 + a 2 )5(p = g e 59 + g v 5y. 


(4.3) 


Here, / and g are the right hand sides (at p = 0) of the first and second equations in (12.5ft . 
respectively, and the indexes 9 and (p denote differentiation with respect to these variables. 
To calculate Ai, we first solve Eqs. (12.51) on the interval (0, t s i m ). Then, using the same initial 
conditions for 59 and 5cp at i — t n (t n — (n — 1)Ar, n = 1, k), we solve Eqs. (14.31) on the 
intervals (t n ,t n + Ar) and from Eq. (14.21) End Ai. 

The reduced energy loss and the lines on which the largest Lyapunov exponent Ai changes 
sign, obtained for At = 2 ■ 10 -5 , Ar = 2 • 10 -3 , k — 1.5 • 10 7 , 59\i n = 1 and 5(p\f n = 0 (this 
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is the usual choice of the initial conditions for 66 and Sip), are shown in Fig. [TJ We note 
three important features of these results. First, the regions in the h-u plane with regular 
(Ai < 0) and chaotic (Ai > 0) dynamics of m are distributed very unevenly. We remark in 
this context that if a contour connecting two points on this plane crosses the white lines an 
even (odd) number of times, then the character of the magnetic moment dynamics in these 
points is the same (different). It should be emphasized, however, that the results concerning 
the regions with regular and chaotic dynamics of the magnetic moment should be considered 
as preliminary. The reason is that, due to the existence of two equilibrium directions of the 
magnetic moment in uniaxial nanoparticles and rather large 86 , the condition Ai > 0 may 
be expected to hold for some regular trajectories as well. In other words, the condition 
Ai > 0 may appear as an artifact of the numerical scheme. The analysis of the long-time 
behavior of m confirms the existence of such trajectories in some regions of the h-u plane 
(these regions are not shown in Fig. [71). Second, the reduced power loss can experience an 
abrupt change not only under transitions between regular and chaotic regimes of precession, 
but also under transitions between different regimes of regular precession (see also Fig. [2j) 
and between different regimes of chaotic precession. And third, comparing in Figs. [3] and 
0 the values of the reduced power loss, one can conclude that nanoparticle heating in the 
circularly polarized magnetic held, whose amplitude and frequency are close to the line of 
P-Q transitions, is more efficient than in the linearly polarized held. On the other hand, 
if the reduced amplitude h is large enough, then the linear polarization of the alternating 
magnetic held is more preferable for heating purpose. 


V. SUMMARY AND CONCLUSIONS 

Using the deterministic Landau-Lifshitz-Gilbert equation, we have studied in detail the 
dependence of the reduced power loss and the character of forced precession of the nanopar¬ 
ticle magnetic moment on amplitude and frequency of circularly and linearly polarized mag¬ 
netic helds. The circularly polarized held, whose plane of polarization is perpendicular to 
the anisotropy axis, can generate three periodic regimes of forced precession (two of them 
occur in the up state of the magnetic moment, and the other occurs in the down state) 
and one quasiperiodic regime. We have determined the regions in the amplitude-frequency 
plane where these regimes exist and have calculated the power loss inside them. A remark- 
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FIG. 7. (Color online) Color map for the reduced power loss Q and the regions of regular (1) and 
chaotic (2) precession of the magnetic moment driven by the linearly polarized magnetic field. The 
change of sign of the largest Lyapunov exponent is indicated by white lines. The parameters p, H 
and a are the same as in Fig. O 

able feature of the power loss is that it changes abruptly at some boundaries between these 
regions. In particular, the transition from the regions with periodic or quasiperiodic preces¬ 
sion to the region with periodic precession in the down state is accompanied by an abrupt 
decrease of the power loss. In contrast, the transition between regions with different periodic 
precessions in the np state is accompanied by an abrupt increase of the power loss if the 
precession angle increases under transition. We have also established that the power loss 
reaches the largest values near the boundary between the region with periodic precession of 
the magnetic moment in the up state and the region with quasiperiodic precession. This 
is the condition under which the nanoparticle heating by the circularly polarized magnetic 
field is the most efficient. 

The linearly polarized field, whose axis of polarization is perpendicular to the anisotropy 
axis, can induce both regular and chaotic regimes of precession of the magnetic moment. 
By analyzing the largest Lyapunov exponent and the long-time behavior of the magnetic 
moment, we have delimited the regions in the amplitude-frequency plane where the magnetic 
moment exhibits the regular and chaotic behavior. The distribution of these regions has a 
complex character and the power loss corresponds, in general, this distribution. Nevertheless, 
the transitions between different regimes of regular and chaotic precession can also strongly 
affect the power loss. Thus, our results provide evidence that the energy dissipation in 
single-domain nanoparticles crucially depends on the character of the magnetic moment 
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dynamics. 
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